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Abstract 

h> ■ 

In general-mass variable-flavor-number (GM-VFN) scheme, the 'heavy quark mechanisms' of 
q ■ B C {B*) meson hadroproduction via the sub-processes g + c — > B C (B*) + • • • and g + b — > B C (B*) H 

on. 

are investigated. In the scheme, possible double counting from these mechanisms and the 'light 

in 

mechanisms' (the gluon-gluon fusion and quark-antiquark annihilation mechanisms) is subtracted 
properly. The numerical results show that the transverse momentum, pt, distribution of the 
Qh. produced B C (B*) from the 'intrinsic heavy quark mechanisms' (i.e. the heavy quark mechanisms in 

which the double counting components have been subtracted accordingly) declines faster with the 
increment of the pt than that from the 'light mechanisms', and only in small pt region (pt ^ 7.0 
GeV) the 'intrinsic heavy quark mechanisms' themselves may make remarkable contributions. The 
combined contributions from the 'intrinsic heavy quark mechanisms' and the 'light mechanisms' to 
the production are compared with that obtained by the most calculations in literature, which in 
some sense are within the fixed flavor number (FFN) scheme at leading order, and we find that the 
production by virtue of the GM-VFN scheme is more or less the same as the one in literature, except 
in the small pt region. 
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I. INTRODUCTION 



The unique and stable 'double heavy-flavored' meson, B C} has been confirmed in experiment 

QBQ 

, and the observations are consistent with the theoretical expectations within theoretical 
uncertainties and experimental errors. In view of the prospects of B c production at Fermilab 
Tevatron (Run-II) and at LHC, the B c physics is compelling. The future copious data require 
more accurate theoretical predictions, especially, that on the production (at Tevatron^_and 
LHC). The hadronic production of B c meson has been studied quite a lot 
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BCVEGPY 
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14 already. It is remarkable that with special techniques the event generator 
12 1 for the production is in compliment to the PYTHIA environment 15] and 



powerful enough in generating event samples for most purposes, i.e., with it one can enhance 
the event generating efficiency greatly in contrast to PYTHIA itself. With newly upgraded 
version BCVEGPY2.0, not only the ground state of B c meson but also 'low laying' excited 
states can be generated. 

Up to now, the predictions for the B c hadronic production are mainly based on the dominant 
gluon-gluon fusion mechanism in terms of perturbative QCD (pQCD) calculations in order of 

, i.e., via the sub-process g + g — > (cb) + b + c with (c6) in the configurations of S- or P-wave 
states. The less important mechanism via light quark- ant iquark annihilation q + q — > (cb) + b+c 
is studied in P, only for comparison. However, the approaches to the existent estimates, in 
fact, all are in the fixed flavor number (FFN) scheme only with certain extension 1 . 

Recently, as pointed out by the authors of Ref. in hadronic production of charmonium, 
the contribution from the sea charm quark via the subprocess c + g — > (cc) + c may be greater 
than that from the gluon-gluon fusion via the subprocess g + g — > (cc) + c + c for charmonium 
color-singlet production. Therefore, it is interesting to examine t he p roduction by applying 
the general-mass variable-flavor-number (GM-VFN) scheme 20j, |2lJ, |22| and see how important 
are the heavy charm and bottom quark mechanisms in the hadronic production of B C (B*) 
meson not only to compare with the gluon-gluon fusion mechanism but also with the existent 
estimations. To be noted that here the heavy quar k means the sea parton for the heavy flavor, 
which is different from what it means in Ref. |26(. Generally the mechanisms induced by the 



1 As for the exact FFN scheme for ^production, the active quark flavor numbers in the initial state is limited 
to Nf = 3 (i.e. the initial state does not contain the heavy quark components at all) and accordingly in the 
production of B C (B*), the PDFs, just CTEQ5F30. 

which corresponds to the 3-flavor scheme of FFN, should 
be used; but in the existent estimates when extending to higher anergy and higher p t for the production, 
CTEQ6L (or CTEQ4L) etc ^25], which is derived with Nf > 3 as the PDFs, was used instead, thus here we 
call such approach as extended FFN scheme. 



heavy quark components must be considered in the GM-VFN scheme, and they are suppressed 
in comparison with the light partons such as the light valance quarks, light sea quarks and 
gluons in PDFs. However, the suppression fact may be 'compensated' by their lower order 
nature in perturbative QCD and a much 'greater' phase space. Naively, the sub-processes 
c + g — > (cb) + b and b + g — > (cb) + c are 2-body — > 2-body processes in the order of 0(a 3 ), 
while the gluon-gluon fusion subprocess g + g — > (cb) + b + c is 2-body — > 3-body process in 
the order of 0(a\). Thus, we devote this work to estimate the production in terms of the GM- 
VFN scheme by taking into account the contributions from the 'heavy quark mechanisms' and 
that from the light partons, especially, the gluon-gluon fusion mechanism together, not only 
estimating the total cross sections, but also studying their properties on transverse-momentum 
p t and rapidity rj distributions etc. 

In GM-VFN scheme when we talk about the heavy quark components of PDFs and taking 
into account both of the 'heavy quark mechanisms' and the gluon-gluon fusion mechanism 
for the hadronic production, one has to solve the double counting problem: i.e. a full QCD 
evolved 'heavy quark' charm/bottom distribution functions, according to the Altarelli-Parisi 
equations, includes all the terms proportional to In (j^r^J (p> the factorization scale and mg the 
heavy quark mass); and some of them come from the gluon-gluon fusion mechanism, i.e., a few 
terms appear from the integration of the phase-space for the gluon-gluon fusion mechanism. 
Therefore, one needs to make proper subtractions to solve the double counting problem j^. 
One convenient way to do the 'subtraction' is to adopt the GM-VFN scheme, in which the 
heavy-quark mass effects are treated in a consistent way both for the hard scattering amplitude 
and the PDFs 21, 3- Moreover, it will be necessary to use the dedicated PDFs with heavy- 
mass effects included, which are determined by global fitting utilizin g m assive hard-scattering 
cross-sections. For instance, for the present analysis, the CTEQ6HQ [33| is used. Later on for 
convenience, we will call the 'heavy quark mechanisms' which have been subtracted according 
to method in GM-VFN scheme as 'intrinsic ones' accordingly. 

The paper is organized as follows. In Sec. II, the basic formulae for the 'intrinsic charm 
and bottom mechanisms' in the GM-VFN scheme are presented for the hadronic production of 
S- wave cfr-quarkonium states, i.e. B c ( 1 Sq) and B*( 3 Si). In Sec. Ill, we present the numerical 
results for the 'intrinsic mechanisms' and make some comparisons about them. The final section 
is reserved for discussion and summary. 



II. CALCULATION TECHNIQUE 



To study the 'intrinsic heavy quark mechanisms' for the hadronic B C (B* )-production, we 
need to study the hadronic processes of g + c — ► B c + ■ ■ ■, g + b — > £> c + ■ ■ • and c + b — >• £? c + ■ • -. 
And, to make the 'intrinsic mechanisms' and the gluon-gluon fusion mechanism consistently 
coexist, we need to do the subtractions so as to avoid the double counting. For this purpose 



we adopt the proper way in GM-VFN scheme 
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in this work. 
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the cross- 



According to pQCD factorization theorem, in the GM-VFN scheme 
section for the hadronic production of B C (B*) (including all of the mechanisms, even the light- 
quark annihilation one via the subprocesses + % — > B c + c + b (k = u, d, s)) is formulated 

as 

da = Ff Ii (x ll ^)Fj {2 (x2 ) ^)®dagg^ Bc{B , ) (x 1 ,X2 ) v) 

+ E F ^(x 1 ,fi)F^(x 2 ,^)(^da qk g- k ^ BciB » ) (x 1 ,x 2 ,fi) 

+ E F 9 H .{x u n) [^02, m) - i^.(x 2 ,/x)(g)F g (x2,M)] <g)da gc _> Bc(B * ) (x 1 ,x 2 ,n) 
+ E F k ( x i,v) [Fh, [x 2 , fi) - F 9 H . (x 2} fi) ® Fg(x 2 , fi)] (g) da g - b ^ Bc(B;) (x 1: x 2 , fi) 
+ E [(F^Xu^-F^Xu^^F^^)) (Fl 3 (x 2 ^)-F 9 Hj (x 2 ^)<^F^x 2 ^))' 

®^e*-»B c (B*)(^l,^2,/i)H , (1) 

where the ellipsis means the terms in higher a s order. Fjj(x,/i) (with H = H\ or #2 and 
x = Xi or x 2 ) is the distribution function of parton i in hadron H. da stands for the hadronic 
cross-section and da stands for the corresponding subprocesses. For convenience, we have taken 
the renormalization scale /j, b for the subprocess and the factorization scale /xp for factorizing 
the PDFs and the hard subprocess to be the same, i.e. /jr = fip = /i. In the square brackets, 
the subtraction term for F®(x,[i) is defined as 

F$(x,iA) SUB = F%{x,ri®Ff{x,ri= f F}{K^)F a H {-^)-. (2) 

y Jx y K K 

The quark distribution F®(x, fi) (with Q stands for heavy quark c or b) within an on-shell gluon 
up to order a s is connected to the familiar g — > splitting function P g ^Q, i.e. 



J?(x jA i) 



a s (/x) /i 



2tt 



In 



/??, 



_ p 

2 ^9 



(3) 



with P„ 



[x 



Hi -2x + 2x 2 ). 



In Eq.fllJ, the first term represents for the dominant gluon-gluon fusion mechanism; the 
second one for the light quark and anti-quark annihilation mechanism; the remainders for the 



FIG. 1: Typical Feynman diagrams for bottom induced process: g{p2) +b(px) — » B c {p^) + c{p^). The 
Feynman diagrams for charm induced process: g + c — > B c + b can be obtained by the replacements: 
b — ► c and c — > 6. 
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FIG. 2: Graphical representation for the subtraction method within GM-VFN I2QH. The subtraction 
term (diagram) is placed in the middle to emphasize its similarity both to the bottom mechanism 
(left) and the gluon-gluon fusion mechanism (right). The symbol x on the internal quark line in 
the subtraction term indicates it is close to the mass-shell and collinear to the gluon and hadron 
momentum. (As 'named' in the text, the first term and the second one with minus sign are considered 
as a whole and called as 'intrinsic one'.) 




12| may be employed 



so 'intrinsic charm/bottom mechanisms', in which the subtraction is introduced to avoid the 
|2fjj . The gluon-gluon fusion mechanism and the light juark and^ anti-quark 

HQ, 



double counting 

annihilation mechanism have been considered in several previous papers 
and for their investigation, now the B c meson generator BCVEGPY [11 , 
conveniently. As for the 'intrinsic charm/bottom mechanisms', we need to calculate three kinds 
of subprocesses: g + c — > B C (B*) +b, g + b — > B C (B*) + c and c + b — > B C (B*). For the hadronic 
production via c + b — > B C (B*), i.e. the fifth term in Eq.(^) when a p t cut (not too tiny) is 
put on (as usually done in experiment), it makes no contributions to the B C (B*) production, 
because its hard subprocess is a 2 — > 1 process essentially. Thus, we will not consider it in this 
paper, because the Tevatron and LHC experiments always put on some cut for small pt- That 
is, below we will concentrated merely on the 'intrinsic heavy quark mechanisms' of B C (B*) 
production via the two sub-processes: g + c — > B C (B*) + b and g + b —> B C (B*) + c. 

Typical Feynman diagrams for these two sub-processes at leading order (LO) are shown in 
Fig. ( | T |) . The corresponding diagrams for B* production are similar, and can be obtained by 



directly replacing B c with B* in Fig.(JT}. To be specific, according to Eq. (JTJ) , the inclusive B c 
hadronic production via 'intrinsic charm/bottom mechanisms' can be formulated as, 

da = J2 J " dx \ J dx 2 F t Hi (x 1 , fx F ) x F^(x 2 , IJ)da ij ^ Bc{B * ) x{x u x 2 , p) , (4) 

where i 7^ j and i,j — g,c for 'intrinsic' charm mechanism and i, j = g,b for 'intrinsic' bottom 
mechanism. Here, the heavy quark PDF F®(x,fi) (x = x\ or x 2 , Q = c or b, H = Hi or H 2 ), 
includes the subtraction term F B (x, p)sub as is defined in Eq.® in order to avoid the double 



counting. Using subtraction method to avoid the doub 



and then as the GM-VFN scheme was developed in 
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e counting was pointed out in Ref. 



221 ] . In Fig. P|) . we take the 'intrinsic' 



bottom process (Fig.(^,)) as an example to illustrate this approach graphically, which is similar 



to the case of hadronic production of heavy quarks [2l|. The symbol x on the internal quark 
line in the subtraction terms mean that the heavy quark four momentum squared is on mass- 
shell and moving longitudinally, which is a good approximation when the quark is collinear to 
the gluon, and results in a factor of order a s distribution of a quark in a gluon, like in Eq.®. 

In Eq.(|IJ), daij^B c x(xi,x 2 , fiF, ^) = d&ij-+ Bc x(xi, x 2 , /i) stands for the usual 2-to-2 differen- 
tial cross-section, 

da ij ^ Bc(B . c)x (x 1 ,x 2 , f x)= 2 ll / 2 - W2F ,/ LP»-Pi-P2 • (5) 

Here, pi, p 2 are the corresponding momenta of the initial two partons and p$, p^ are the 
momenta of the final ones, respectively. The initial-parton spin and color average and the 
final-state quantum number summation are all attributed to the |M| 2 . According to GM- 
VFN scheme, the heavy quark masses are kept in the evaluation of S-matrix. For shortening 
the text, we put the explicit expression of amplitude squared, |M| 2 , in the Appendix instead. 
We generate the amplitude and square it |M| 2 with the program package: Feynman Diagram 
Calculation (FDC) 28], which is a Reduce and Fortran package to perform Feynman diagram 



calculation automatically and has been well-tested by various applications. The phase space 
integration is mani pul ated by the routines RAMBOS 29] and VEGAS jsO], which can be found 
in BCVEGPY 
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III. NUMERICAL RESULTS AND DISCUSSIONS 



In this section, we present the numerical results of the 'intrinsic charm and bottom mech- 
anisms' and then make a comparison with that of the gluon-gluon fusion mechanism. The 
combined results for both the 'intrinsic mechanisms' and the gluon-gluon fusion mechanism are 



consistently treated within the GM-VFN scheme. Finally, we shall make a comparison of the 
present results with the existing results in the literature. 



A. The results for the 'intrinsic' charm and bottom mechanisms 

In doing the numerical calculations, the values of the radial wave function of B c or B* at the 



origin is taken as what in Refs. 



31 



, 32|, i.e.. | /?(()) | 2 = 1 ■ ■:> IGV \ the masses of <• and i> <|uarks 
are taken as m c = 1.50 GeV and m& = 4.80 GeV, respectively; the mass of the bound state is 
approximately taken to be the sum of the two heavy quark masses, i.e. M = mi, + m c = 6.30 



GeV (to ensure the gauge invariance of the concerned amplitudes). 

There are a couple of uncertainties remaining in the theoretical estimations for the B c meson 
hadronic production ^(J, such as those from the choice of renormalization scale, factorization 
scale, etc.. In present calculations, the factorization and the renormalization energy scales are 



set to be the 'transverse mass' of the bound state, i.e. Q = M± = y M 2 + p 2 , where p t is 
the transverse momentum of the bound state; here CTEQ6HQ [3J| for PDF and the leading 
order a s running above ^qcd^ = 0.326 GeV are adopted. CTEQ6HQ is adopted here in 
the calculations, because it is an improved set of the parton distributions determined for the 
GM-VFN scheme that incorporates heavy flavor mass effects. 

To see the fact how the 'double counting' is subtracted in GM-VFN scheme precisely, first 
of all we take LHC as an example to compute each term numerically according to Eq.(pQ) for 
heavy quarks and present them in Fig.©. From the figure, one may see the consequence of 
subtraction of the 'double counting' quite substantially. There is large cancellation between the 
contributions from the 'pure' terms (with heavy quark's PDF taken to be CTEQ6HQ) and the 
subtraction terms (with heavy quark's PDF taken to be the subtraction terms defined in Eq.(J2J)), 
especially, in the large p t regions. Below when we compare the GM-VFN scheme with FFN 
scheme for the production, we will, in fact, compare the results for the production obtained 
by sum of the gluon-gluon fusion mechanism and the 'intrinsic heavy quark mechanisms' in 
GM-VFN scheme with those obtained by gluon-gluon fusion mechanism with PDFs extended 
to CTEQ6L in FFN scheme. 

In high energy hadronic collisions, events with a small pt and/or a large rapidity y are hard to 
be detected by detectors, or say very difficult to be reconstructed among the backgrounds. So in 
experimental observations, the events with small pt and large rapidity y are dropped practically. 
Therefore, in theoretical estimates, proper cuts onp t and y are applied. TABLEs HI and HP show 
the cross-sections for the hadronic production of the B C (B*) meson at LHC and TEVATRON 
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FIG. 3: The p^-distributions of the 'intrinsic mechanisms' of the hadronic production of B c at the 
LHC. The left diagram is for g + b — > B c + c and the right diagram is for g + c — > B c + b, where the 
'PURE-PDF' ones mean the heavy quarks' PDFs are taken just to be CTEQ6HQ, 'SUB' ones mean the 
heavy quarks' PDFs are taken just to be the subtraction term defined in Eq.Q, while 'COMBINED' 
ones are the combination of the 'PURE-PDF' and 'SUB' components with proper sign as indicated by 
Eq.([T|>. which correspond to the 'intrinsic charm g + c — > B c + b or bottom g + b — > B c + c' precisely. 

TABLE I: The cross-section (in unit of nb) for the hadronic production of B c at LHC (14.0 TeV) and 
TEVATRON (1.96 TeV), where for simplicity the symbol g + b means g + b — > B c + c and etc. In the 
calculations, pt > AGeV is taken. \y\ < 1.5 for LHC, while \y\ < 0.6 at TEVATRON. q stands for the 
sum of all the light quarks (u, d and s). 
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3.28 x 10~ 3 


0.542 
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respectively with cuts: p t > 4GeV for both LHC and TEVATRON, and \y\ < 1.5 for LHC, 
\y\ < 0.6 for TEVATRON. In TABLEs |U and |Hj we also show the results from CTEQ6L 
which will be discussed later. The tables show that the cross-sections of the 'intrinsic charm and 
bottom mechanisms' are comparable to the gluon-gluon fusion mechanism, and the sequential 
order for the cross-sections is a gg > cr g i > a gc » a q q. Because the contributions from the 
light quark and anti-quark annihilation mechanism are quite small comparing to the other 



TABLE II: The cross-section (in unit of nb) for the hadronic production of B* at LHC (14.0 TeV) 
and TEVATRON (1.96 TeV), where for simplicity the symbol g + b means g + b — > B* + c and etc. In 
the calculations, pt > AGeV is taken. \y\ < 1.5 for LHC, while \y\ < 0.6 at TEVATRON. q stands for 
the sum of all the light quarks (u, d and s). 
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FIG. 4: The pt distributions of the hadronic production of B c (left diagram) and B* (right diagram) 
at LHC in GM-VFN scheme. The solid line is for the g + g —> B C (B*) mechanism, the dash line is 
for the intrinsic g + b — > B C (B*) mechanism and the dotted line is for the intrinsic g + c — > B C (B*) 
mechanism. All the pt distributions are drawn with \y\ < 1.5 and the PDF is taken to be CTEQ6HQ. 

mechanisms, below we will neglect them. 

The contributions from the 'intrinsic charm and bottom mechanisms' are remarkable in small 
Pt region of the B C (B*) production. To show this point clearly, we present the the transverse 
momentum distributions of B C (B2) in Figs. (JU EJ) . From Figs. 003), one may observe that in the 
small pt regions, the 'intrinsic charm and bottom mechanisms' are comparable or even greater 
than that of the gluon-gluon fusion mechanism. One point worthy to mention is that even 
though the 'intrinsic charm mechanism' are smaller than that of 'intrinsic bottom mechanism' 
for the production in the small pt region for both B c and B* production, whereas in the high p t 




p t (GeV) 



p t (GeV) 



FIG. 5: The pt distributions for the hadronic production of B c (left diagram) and B* (right diagram) 
at TEVATRON in GM-VFN scheme. The solid line is for the g+g — > B C (B*) mechanism, the dash line 
is for the intrinsic g + b — ► B C (B*) mechanism and the dotted line is for the intrinsic g + c — > B C (B*) 
mechanism. All the pt distributions are drawn with \y\ < 0.6 and the PDF is taken to be CTEQ6HQ. 

region, the 'intrinsic charm' contribution becomes slightly larger than that of 'intrinsic bottom' 
contribution for the case of B*. The contributions from both 'intrinsic charm and bottom 
mechanisms' drop much more rapidly with the increment of pt than those from the gluon-gluon 
fusion mechanism. 



B. Comparisons between the results 

In the literature, only the gluon-gluon fusion mechanism for the hadronic production of 
B C (B*) is considered in most of the calculations, while the light quark and anti-quark anni- 
hilation mechanism, being less important, is dropped. In these calculations, which can be 
considered as an extension of the FFN scheme as explained at the above, therefore there is 
no 'intrinsic charm/bottom mechanisms'. In the FFN scheme, the active partons in the initial 
state are limited: only light quarks rif = 3 and gluons, while the heavy charm and/or bottom 
quarks appear only in the final state. To be consistent with the exact FFN scheme, the PDFs 



for the initial partons should be CTEQ5F3 [241 . the one generated by using the evo 



ution 



ker- 



nels with effective flavor number n e // = 3. However, in the references 
to do t 
in Ref. 

very great (< 10%). In fact, our numerically calculation shows that it only leads to a small 
difference ~ 5% with CTEQ6L to replace CTEQ5F3 in the estimates. 



le production estimates, CTEQ6L for PDFs is used instead of CTEQ5F3. It is shown 
Tol | that generally the uncertainties caused by different leading order PDFs cannot be 




FIG. 6: The pt distributions of the hadronic production of B c (left diagram) and B* (right diagram) at 
LHC. The solid and the dashed lines stand for the total (sum of the 'intrinsic heavy quark mechanisms' 
and the gluon-gluon fusion mechanism) results obtained by the GM-VFN scheme for rapidity cuts 
\y\ < 2.5 and \y\ < 1.5 respectively; and the dash-dot and the dotted lines for gluon-gluon fusion 
results obtained by the FFN scheme with PDFs CTEQ6L for rapidity cuts \y\ < 2.5 and \y\ < 1.5 
respectively. 

In TABLEs U] and |H] the cross-sections obtained with GM-VFN scheme and FFN scheme 
(with slight extension CTEQ6L to replace CTEQ5F3) for the hadronic production of B C (B*) 
are shown quantitatively. Due to the fact that the gluon distribution of CTEQ6HQ for GM- 
VFN scheme is always smaller than that of CTEQ6L, especially at small x regions, so the cross- 
section for the gluon-gluon fusion under the GM-VFN scheme is smaller than that under the 
FFN scheme. Moreover since x may reach to much smaller region at LHC than at TEVATRON 
for the production, the difference between these two schemes is bigger at LHC than that at 
TEVATRON. TABLEs HI ITU shows this point clearly: at LHC, the cross-section for gluon-gluon 
fusion under the GM-VFN is only ~ 60% of the case of FFN, while at TEVATRON, such ratio 
is changed to ~ 80%. When taking the 'intrinsic heavy quark mechanisms' into account for 
the GM-VFN scheme, one may find that the gap between the GM-VFN results and the FFN 
results can be shrunk sizably. 

In fact, the difference of the predictions by the GM-VFN scheme and the extended FFN 
scheme is mainly in small p t region. When p t becomes 'big', according to the separation 
of the GM-VFN scheme, the contributions from the 'intrinsic heavy quark mechanisms' are 
greatly suppressed due to the cancellation from the subtraction terms so as to avoid the double 
counting with those from the gluon-gluon fusion mechanism, hence the contributions from the 
gluon-gluon fusion mechanism become dominant in the GM-VFN scheme. The sum of the con- 
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FIG. 7: The pt distributions of the hadronic production of B c (left diagram) and B* (right diagram) 
at TEVATRON. The solid and the dashed lines stand for the total (sum of the 'intrinsic heavy quark 
mechanisms' and the gluon-gluon fusion mechanism) results obtained by the GM-VFN scheme for 
rapidity cuts \y\ < 1.5 and \y\ < 0.6 respectively; and the dash-dot and the dotted lines for gluon- 
gluon fusion results obtained by the FFN scheme with PDFs CTEQ6L for rapidity cuts \y\ < 1.5 and 
\y\ < 0.6 respectively. 

tributions from the 'intrinsic heavy quark mechanisms' and the gluon-gluon fusion mechanism 
in the GM-VFN scheme happens to close to those from FFN scheme, where only the gluon- 
gluon fusion mechanism alone is dominant. In order to see the fact precisely, we present the p t 
distributions predicted by the two schemes for the hadronic production of B C (B*) at LHC and 
TEVATRON with two possible rapidity cuts in FIGs.(JBJ) and ((Zj) respectively. From FIGs.© 
and (J7|), one may observe clearly that the main difference between the predictions by GM-VFN 
scheme and the extended FFN scheme (in literatures) is only in small pt region (p t ^ 6.0 ~ 7.0 
GeV). In higher p t regions (pt > 6.0 ~ 7.0 GeV), the curves predicted by the two schemes 
become very close to each other. 

IV. DISCUSSION AND SUMMARY 

We have studied the hadronic production of B C (B*) meson induced by the heavy quarks 
inside the incident hadrons under the GM-VFN scheme. The double counting problem of the 
gluon-gluon mechanism with the 'intrinsic heavy quark mechanisms' is consistently solved in 
the GM-VFN scheme, and also the heavy quark mass effects, which are important for precise 
investigations, are taken into account in both heavy quark PDFs and the hard scattering kernel. 

As shown in TABLEs (jlj) and (JTl]), we have found that the total cross sections of the 'intrin- 



sic charm and the bottom mechanisms' may be comparable to that of the gluon-gluon fusion 
process, especially the bottom quark interaction mechanism. The B c production from the 'in- 
trinsic heavy quarks' drops much more rapidly with the increment of the transverse momentum 
p t than that from the gluon-gluon fusion mechanism, although it may 'overrun' the later in the 
small transverse momentum region. At about 6.0 ~ 7.0 GeV the differential cross section of 
the B c transverse production rate due to the gluon-gluon fusion mechanism surpasses all that 
from the involved 'intrinsic charm and bottom mechanisms'. 

The calculations show clearly that the difference between the GM-VFN scheme and the FFN 
scheme (adopted in most of the calculations in the literature) is mainly around the small pt 
regionp t < 6.0 ~ 7.0 GeV, as shown in Figs.© and (J7J). While the value of p t becomes big, 
the results are almost the same for both schemes. This is because that when the magnitude 
of pt becomes big, the contributions from the 'intrinsic heavy quark mechanisms' are greatly 
suppressed due to the subtraction terms in GM-VFN scheme and the dominant contribution 
tends to be the gluon-gluon fusion mechanism. For the most practical purposes and for the 
production at LHC and Tevatron one does not need to care of the difference on the predictions 
of the GM-VFN scheme and the extended FFN scheme which were used in most references 
except those in the small region (p t < 6 ~ 7 GeV). 

In conclusion, if one takes the GM-VFN scheme, the 'intrinsic heavy quark mechanisms' 
in B C (B*) hadroproduction are important in small p t region in comparison with the dominant 
gluon-gluon fusion mechanism. Therefore, if pt-cut of B c can be taken so small as 4 GeV, the 
contributions from the 'intrinsic heavy quark mechanisms' may be seen by the difference from 
the predictions not only of the gluon-gluon fusion mechanism in GM-VFN scheme, but also of 
the mechanism in extended FFN scheme (the predictions in most references). Hence a detailed 
experimental study of the production in small pt region could gain some information on the 
'intrinsic charm and bottom' distributions inside the hadron, and will give a 'judgement' about 
the GM-VFN scheme and the extended FFN scheme. According to the investigation here, 
the 'intrinsic heavy quark mechanisms' can be used as a supplement to the usual gluon-gluon 
fusion mechanism in GM-VFN scheme, especially, when experiments can reach to very small 
p t region and indicate some deficiency for the production. Probably some suitable fixed target 
experiments, in which the detector may cover almost all solid angles (almost without p t cut), 
can test the 'intrinsic charm and bottom mechanisms' in the future. 
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Appendix: The squared amplitudes for the sub-processes 

For convenience, we express the square of the amplitudes by the Mandelstam variants s, t 
and u, which are defined as: 

s = (pi + P2) 2 , t=(p 1 - p 3 ) 2 , u =(p 1 - p A ) 2 , 

where pi = (Ei,pi X ,piy,p iz ) are the corresponding momenta for the involved particles: p\ and]?2 
are the momenta of initial partons, p 3 and p 4 are the momenta of B C (B*) and another outgoing 
particles respectively. Further more, for b(pi) + g{p2) — > B c (p 3 )/ B*(p 3 ) + c{p A ), we set 

Mi = (u - (m b + m c ) 2 ), si = (s-m 2 b ), t ± = (t-m 2 c ). 

The relation, u\ + t± + si = 0, is useful to make all the expressions for the square of the 
amplitudes compact. 

The square of the amplitude for the subprocess b(pi) + g{p2) — > B c (p 3 ) + c(p A ) can be written 

as, 
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where (s x <-> ti, m fe <-> m c ) stands for the remaining terms that can be directly obtained by 
exchanging si <-> ti and mj <-> m c for all the terms in the first big parenthesis. The square of 
the amplitude for the subprocess c(pi) + g(j>2) — > B c (p 3 ) + 6(1)4) can be directly obtained from 
the above formula by taking the transition <-> m c . 
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where (s% <-> ti, <-> m c ) stands for the remaining terms that can be directly obtained by 
exchanging si «-> ti and <-> m c for all the terms in the first big parenthesis. The square of 
the amplitude for the subprocess c(pi) + g{p2) — > B*(p 3 ) + b{p^) can be directly obtained from 
the above formula by taking the transition m b «-> m c . 
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